Role of heating and current-induced forces in the stabihty of atomic wires 
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We investigate the role of local heating and forces on ions in the stability of current-carrying 
aluminum wires. We find that heating increases with wire length due to a red shift of the frequency 
spectrum. Nevertheless, the local temperature of the wire is relatively low for a wide range of biases 
provided good thermal contact exists between the wire and the bulk electrodes. On the contrary, 
current-induced forces increase substantially as a function of bias and reach bond-breaking values 
at about 1 V. These results suggest that local heating promotes low-bias instabilities if dissipation 
into the bulk electrodes is not efficient, while current-induced forces are mainly responsible for the 
wire break-up at large biases. We compare these results to experimental observations. 



Atomic wires are an ideal testbed to study transport 
properties at the nanoscale P, |3| ■ Physical phenomena 
that have been investigated in these systems include their 
quantized conductance conductance and noise oscil- 
lations as a function of the wire length |3, S Ifil 13 ^heat- 
ing 1^ and current-induced atomic motion 0, 
The latter two properties, in particular, are of 
key importance in understanding the stability of atomic 
wires under current flow |l2l| . 

Halbritter et. al. and, more recently, Mizobata et. 
al. 0, 01 have explored the mechanical stability of 
atom-sized Al wires. These authors found that the proba- 
bility of forming single-atom contacts decreases with in- 
creasing bias and that it vanishes at a critical bias of 
about 1 V In addition, several samples exhib- 

ited low-bias instabilities (typically at biases of 100 mV 
or less) leading to break-up of the atomic wires. Sim- 
ilar trends have been reported in the case of Pb [l^ 
and Au ^ point contacts, suggesting that the mecha- 
nisms leading to such instabilities are material indepen- 
dent. These low- bias instabilities have typically been at- 
tributed to inadequate dissipation of heat into the bulk 
electrodes. 4, 8, 9] However, additional forces on ions due 
to current flow may also contribute, since both effects are 
present at any given bias. 



In this Letter, we explore the relative role of local 
heating and current-induced forces in the stability of alu- 
minum wires at different biases. We study these effects 
perturbatively, i.e., we first calculate current-induced 
forces on ions while neglecting local heating; second, we 
assume current-induced forces are negligible, and evalu- 
ate the local temperature of the wires. This perturbative 
approach is supported a posteriori: we find that local 
heating is substantial at very low biases if poor ther- 
mal contacts exist between the wire and the bulk elec- 
trodes H, 0, 01 , while current-induced forces (in partic- 
ular, the average force per atom, see below) are small 
at low biases but increase substantially with increasing 
external voltage. In addition, if the heat in the wire 
can be dissipated efficiently into the electrodes, the local 
temperature in the junction will be quite low at biases for 
which current-induced forces reach bond-breaking values. 
These results suggest that low-bias instabilities of atomic 
wires, if any, are due to local heating while high-bias 
break-up of a junction is mainly due to current-induced 
forces. 

A schematic of one of the wires investigated is depicted 
in the inset of Fig. ^a). It consists of an Al atomic 
chain (the inset shows a one-atom wire) sandwiched be- 
tween two Al electrodes that we model with ideal met- 
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als (jellium model, rs R:! 2) [UliliJ. The stationary 
scattering wave functions of the whole system are calcu- 
lated by solving the Lippmann-Schwinger equation iter- 
atively to self-consistency. Exchange and correlation are 
included in the density-functional formalism within the 
local-density approximation [19L Eol l2l| . More details 
of the calculations of current, forces and heating can be 
found in the original papers [sl, 113; ■ 

We investigate the transport properties of Al chains 
with up to four atoms. All atomic coordinates have been 
relaxed to the equilibrium value at zero bias |2,ll |. The 
relaxed Al-jellium edge bond length is about 2.0 a.u. and 
the relaxed Al-Al bond distance decreases from 6.8 a.u. 
in the two-Al wire to about 5.8 a.u. in the four-Al wire. 
The bond lengths in longer wires are shorter and ap- 
proach the bulk bond length, due to the reduced interac- 
tion of the middle atoms with the electrodes. 

Current-induced forces - The force on an ion, with po- 
sition R, due to current flow is calculated as [13, 123, l2J| 
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where v is the electron-ion interaction. The first term 
on the RHS of Eq. ^ is similar to the usual Hellmann- 
Feynman contribution to the force due to localized elec- 
tronic states {ipi)- The second term is the contribution 
due to continuum states jl^. The wavefunctions |^a} 
are eigendifferentials for each energy interval A in the 
continuum a j2^ . Finally, the forces due to ion- ion in- 
teractions are added to obtain the total force. The total 
current-induced forces as a function of applied bias in 
Al wires containing one, two, three and four atoms are 
shown in Fig. Q^a), (b), (c), and (d), respectively. The 
atoms in the chains are denoted as Al^ {i — 1,...,4) 
starting from the left electrode. Positive force indicates 
atoms are pushed to the right, i.e., against the electron 
flow. For the small external bias of 0.01 V, the forces on 
All and AI2 in the 2-Al wire are in the same direction 
and have almost equal values (see Fig. H^b)), indicat- 
ing the symmetry of the two atoms. For 3-Al and 4-Al 
wires, the forces obey the zero-sum rule as also found in 
Refs. [lOl, |25f the force on AI2 is opposite to that of Ali 
and AI3 in the 3-Al wire, and the ones on Ali and AI4 
are opposite to those of AI2 and AI3 in the 4-Al wire, and 
the sum of the forces of each wire is close to zero. 

For 1-Al and 2-Al wires current-induced forces move 
the atoms against the electron flow at any bias (see Fig.nj 
and increase quite linearly with increasing voltage up to 
IV. On the other hand, for the 3-Al and 4-Al wire cases, 
the forces on ions are nonlinear at small voltages with 
some atoms moving against and some along with the cur- 
rent flow. The direction and magnitude of these forces is 
related to the local resistivity dipoles that form around 
the atoms |0. 
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FIG. 1: Total current-induced forces as a function of bias in 
atomic wires containing (a) one, (b) two, (c) three, and (d) 
four atoms. The inset shows a schematic of one of the wires 
investigated. The atoms are labeled from the left electrode. 
The left electrode is positively biased. Positive force pushes 
the atom against electron flow. 
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FIG. 2: (a) Average force on ions as a function of bias for 
Al wires of different length, (b) Average current per ion as a 
function of bias for the same wires. 
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Despite these nonlinearities as a function of bias, the 
average force per atom follows a simple trend: it increases 
almost linearly with bias and for a given external voltage 
it is smaller the longer the wire (see Fig. |2Ia)). Simi- 
lar results have been reported for Si atomic chains [lo| . 
This trend correlates well with corresponding current per 
atom as a function of bias (Fig.|2Ib)) and it shows that 
both the force and current per atom saturate with length. 
This is consistent with the fact that (i) the wire has no 
defects and the resistance is mainly determined by the 
contacts, and (ii) the average force due to charge redis- 
tribution reaches a constant value with increasing wire 
length, since the stationary states become evenly dis- 
tributed across the length of wire. 

While it is difficult to predict the exact bias at which 
an atomic chain breaks due to current-induced forces, we 
can estimate it using total-energy calculations. For Al 
atomic wires, Jelfnek et al. have found by first-principles 
calculations that the maximum force the wires can stand 
is about 0.6 eV/A 0]. From Fig. EJa) such force cor- 
responds to a bias of about 0.8 V. This value compares 
very well with the reported experimental result of 0.8 V 
by Mizobata et. al. for single-atom Al contacts 0, . 
It is also evident that, on average and everything else 
being equal, larger voltages are required to break longer 
wires (Fig. [Sfa)). However, the break-up of the wires at 
high bias is likely to nucleate from the bonds of the atom 
with the largest force. 

Local heating - We now explore the effect of heating 
assuming current-induced forces are zero. Heating oc- 
curs when electrons exchange energy with the lattice via 
absorption and emission of vibrational modes. Details of 
the theory of local heating in nanoscale structures can be 
found in Refs. 9, 27. Here we just mention that there are 
two major processes that lead to a given local tempera- 
ture in a nanojunction. One is due to inelastic processes 
that occur in the region of the junction. In this case elec- 
trons incident from the right or left electrode can absorb 
(cooling) or emit (heating) energy because of electron- 
vibration scattering in the junction. The other is due 
to dissipation of energy into the bulk electrodes. Let us 
ffi'st focus on the inelastic scattering contribution assum- 
ing the dissipation into the electrodes is negligible. This 
can be the result of, e.g., weak coupling of vibrational 
modes localized in the junction with the continuum of 
modes of the bulk electrodes 1271. 



Denoting by Wi 



the power absorbed (emitted) 



by electrons incident from the left (right) via a vibra- 
tional mode v, the total thermal power generated in the 
junction can be written as the sum over all vibrational 
modes of four scattering processes 0: 



P = 



E 



w, 



R,2 



W, 
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W, 
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This power can be expressed in terms of the electron- 
vibrational coupling in the presence of current. This 
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FIG. 3: Local temperature as a function of bias for Al wires of 
different lengths. No heat dissipation into the bulk electrodes 
is taken into account. 
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FIG. 4: Local temperature as a function of bias for Al wires 
of different lengths. Heat dissipation into the bulk electrodes 
is taken into account. 



coupling, in turn, can be determined from the same sta- 
tionary scattering wave functions used to calculate forces 
on the ions 0. When the heating processes (W^^'^ and 
VF^'^) balance the cooling processes (W^'^ and W^^^), 
i.e., P = 0, a steady-state local temperature is estab- 
lished in the junction. This temperature is plotted in 
Fig. 131 assuming zero background temperature. 

A bias larger than the first vibrational frequency of the 
junction is necessary to generate heat (see Fig. ^ [2^ . 
One can see that the equilibrium temperature increases 
abruptly above that threshold bias and becomes substan- 
tial at biases of only few mV (Fig.OJ. This implies that 
if the heat generated in the junction is not efficiently 
dissipated into the electrodes, very large local tempera- 
tures will be generated even at low biases; thus induc- 
ing structural instabilities and consequent wire break- 
up. These instabilities have indeed been observed in Al 
wires 0,0,0], as well as in other metalHc wires |3.l2^. 
Due to a red shift of the frequency spectrum with increase 
of wire length at a given bias, longer wires have substan- 
tially larger local temperatures. Note that at these same 
biases, current-induced forces are neghgible (see Fig.^). 
However, if there is sufficient coupling between the junc- 
tion vibrational modes and the modes of the bulk elec- 
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trodes, this will allow the heat to dissipate from the junc- 
tion into the bulk electrodes via elastic phonon scatter- 
ing Is^l • This results in a much lower local tempera- 
ture of the junction (see Fig. It is clear that at the 
biases where current-induced forces reach bond-breaking 
values, the temperature in the junction is still quite low 
(less than about 50 K for all wires at 1 V) with shorter 
wires having lower temperatures. 

In conclusion, we find that both local heating and 
current-induced forces can lead to structural instabil- 
ities, but at different ranges of the external bias. If 
there is poor thermal contact between the junction and 
the electrodes, local heating induces instabilities in the 
structure at small biases. At high-biases, break-up of 
nanoscale junctions is mainly due to current-induced 
forces. These results help illuminate a very important 
aspect of nanoscale electronics. 
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deeNL (e) Nr (e) [ul (e) - nR{e)] , (3) 



where nj^j^j) is the Bose-Einstein distribution function and 
Al{_r) (e) is the spectral density of states at the left (right) 
electrode surface. 



